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ABSTRACT 

Nucleosynthesis of heavy nuclei in metal-poor stars is generally thought to 
occur via the r-process because the r-process is a primary process that would have 
operated early in the Galaxy's history. This idea is strongly supported by the fact 
that the abundance pattern in many metal-poor stars matches well the inferred 
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Solar r-process abundance pattern in the mass range between the second and 
third r-process abundance peaks. Nevertheless, a significant number of metal- 
poor stars do not share this standard r-process template. In this Letter we suggest 
that the nuclides observed in many of these stars are produced by the r-process, 
but that it is prevented from running to completion in more massive stars by 
collapse to black holes before the r-process is completed, creating a "truncated 
r-process," or "tr-process." We find that the observed fraction of tr-process stars 
is qualitatively what one would expect from the initial mass function, and that 
an apparent sharp truncation observed at around mass 160 could result from a 
combination of collapses to black holes and the difficulty of observing the higher 
mass rare earths. We test the tr-process hypothesis with r-process calculations 
that are terminated before all r-process trajectories have been ejected. We find 
qualitative agreement between observation and theory when black hole collapse 
and observational realities are taken into account. 

Subject headings: stars: Population II — nucleosynthesis — black hole physics 



Introduction 



The r-process has been understood for many years to (1) synthesize half the nuclides 
heavier than iron, (2) synthesize all of the nuclides heavier than 209 Bi, and (3) be primary 
i n the sense that its nucleosynthesis does not appear to depend on preexisting nuc lides 
feurbidee et~aD Il957l : Wooslev et al1ll994J : IWallerstein et al.lll997t iFarouqi et all I2OO9I ) . In 
this latter context, its productio n of the heavier r-proce ss nuclides in many metal poor 
stars appears to be very uniform ( iSneden fc Cowanl 120031 ) . and it also produces relative r- 
process abundances that are essentially the same as those found in more modern stars. A 
standard interpretation of the r-proces s suggests that it occurs in the neutrino wind that em- 



anates from core-co llapse supernovae (jWoosley et al.lll994j ; iTakahashi. Witti. fc Ja nka 



1996; 



Farouai et al.ll20091). although that interp r etation is not without its is su es (IQian fc Woosley 



1994 



Mever. McLaughlin, fc FuUej Il998t iHudepohl. L. et all hoiOal lbl: iFischer et all boiol ; 



Roberts. Wooslev. k Hoffman! hoioh . Even with these difficulties, thou gh, there are a suf- 
ficient number of uncertainties, e.g., the possibility of sterile neutrinos (IMcLaughlin et al. 



19991 ; IMcLaughlin et al.l 120031 ). that this r-process site might turn out to be viable. 



At the same time, other metal poor stars exhibit some of the features of r-process nu- 
cleosynthesis, but thei r abundances rep r esent a surprisingly poor match with the " standard" 
abundance template ( Aoki et al. 2000 ; Honda et al. 2007 ; Roederer et al. boiOa), p erhaps 
most frequently identified as that observed in CS-22892-052 JSneden fc Cowanl 120031 ). A re- 
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cent paper (IRoederer et al.ll2010al ) has summarized the situation that exists for metal-poor 
stars. The data in that paper seem to suggest that a distribution of r-process abundances 
exists in metal poor stars, with some resembling the standard abundance set, but with a 
significant fraction of stars having abundances that do not match the standard template. 
These latter stars appear to favor the lighter r-process nuclides at varying levels, and many 
seem to have abundance patterns that terminate around Dy, that is, around mass 160. 

In this Letter we point out that the abundance patterns observed for the stars that 
do not fit the standard r-process template could be produced by stars that are sufficiently 
massive that their core collapse first produce neutron stars but that the infall that occurs 
following the formation of the neutron star subsequently results in collapse to a black hole - 
so called "fall back supernovae." This class of stars spans a mass range from roughly 25 to 40 
solar masses (jlieger et al. 1992 ) for low-metallicity stars. When the neutron star collapses 



to a black hole the r-process enrichment of the interstellar medium would cease, terminating 
either when the r-processed regions were s wallowed by the black hole, or when the electron 



antineutrinos fell below the event horizon ( iSasaqui. Kajino. fc Balantekinl 120051 ) . Thus, this 



truncated r-process, or tr-process, nucleosynthesis would terminate at different stages of that 
process, depending on the precise time at which the black hole prevented further r-process 
production or emission of nuclides into the interstellar medium. We suggest that the delayed 
collapse to the black hole, combined with another effect, namely, the difficulties in observing 
the higher mass rare earths, could produce the cutoff in the r-process distributions observed 
around 160 u. 

The general scenario we envision merely assumes that mass elements that predominantly 
produce the lighter r-process nuclei are ejected before the mass elements that produce the 
heavier r-process nuclides. Any setting within a core-collapse supernova that satisfies this 
condition would allow for a truncated r-process. To evaluate this idea quantitatively, how- 
ever, we apply the neutrino-driven wind model for the r-process. Although, as noted above, 
this model is not without issues; it is plausible and well discussed in the literature and makes 
a good setting for discussing the truncated r-process. 



2. Neutrino Wind Model of the r-Process 

In the neutrino wind scenario of the r-process, neutrinos from the nascent neutron star 
heat material at the neutron star surface and drive it away in a wind. As a wind element 
leaves the neutron star surface, it expands and cools. Nucleons in the wind ele ment assemble 



into heavier nuclei, which serve as the seeds for the subsequent r-process (jWoosley et al. 



19941 ) . Importantly, the abundance pattern for the seed nuclei is well- described by a quasi- 



-4- 



equilibrium, which, for the entropy and neutron richness characteristic of neutrino-driven 
wind envirnoments, has a relatively sharp abundance peak near A^IOO. As the temperature 
in the wind element falls further, charged-particle reactions freeze out. If a high neutron 
density accompanies the seed nuclei, neutron captures and beta decays successively promote 
the seed nuclei to higher mass. This upwards flow in mass slows at the neutron closed shells 
at 82 and t hen 126 neutrons, and thereby produces th e r-process abundance peaks at A~130 
and 195 u JWooslev et al.lll994j ; IWanaio et al~ll2002h . 



For what follows, it is essential to note that the neutrino-driven wind leaves the star 



over a period of several seconds. In the standard scenario (IMeyer et al.lll992l ; IWoosley et al. 



19941 ; iTakahashi. Witti. &: Jankal Il994[ ). the wind elements that leave the star early have 
lower entropy and a lower degree of neutron richness. As the neutron star evolves, the wind 
elements that leave later have higher entropy and a greater degree of neutron richness. It 
is in these latter wind elements that the heaviest nuclei are made. The lower-entropy, less 
neutron rich wind elements that leave the neutron star earlier produce ligher nuclei. If the 
neutron star collapses to a black hole after the lower-entropy elements leave the star but 
before the higher-entropy do, then the r-process has been truncated, and the abundance 
pattern will be dominated by lighter r-process nuclei. 

Current supernova models do not naturally produce the parameters required to produce 
a successful r-process, most notably the entropy, but also the requisite neutron richness. 
While it may in fact be the case that core-collapse supernovae are not the site of the r- 
process, it may also be that current supernova models are not sufficiently detailed to provide 
an adequate description of the r-process conditions in core-collapse events. For example, 
results from multidimensional hydrodynamics calculations suggest that the instabilities re- 
sulting from those calculations may ulti mately be shown to p roduce the entropy required 
for making the heaviest r-process nuclei (IBurrows et al.l 120061). Interestingly, these insta- 
biliti es may generate the neutron star kicks ( Burrows et al. 20061; Guilet. Sato, h Foglizzo 



20101 ) that give rise to high space velocities of some pulsars (lArzoumanian et al. 



2002). It 



will also be necessary to include all of the detailed neutrino physics (IBurrows fc Thomson 
20021 ; iRoberts. Woosley. & Hoffman! I201CH : iLiebendoerfer et al.l 120051 ) to adequately charac- 
terize the neutron richness in matter expelled from the supernova. 

Future work will determine whether more advanced core-collapse supernova models will 
provide the conditions needed for the r-process. For the purposes of the present paper, we 
will simply assume that neutrino-driven winds in core-collapse supernovae are at least one 
of the sites of r-process nucleosynthesis. 
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Model Calculations 



To study the trunc ated r-process i n the neutrino-driven wind, we applied the basic 
idea behind the study of IWoosley et al.l (119941 ). which assumed that the r-process occurred 
in the neutrino wind from a core collapse supernova. In that study, a succession of 40 
"trajectories" (that is, thin shell wind elements), all assumed to have originated deep within 
the (assumed spherically symmetric) star, but having initially different conditions of density, 
temperature, entropy, and electron fraction, were emitted into the interstellar medium, thus 
contributing to the total r-process nucleosynthesis. The bubble was evolving in time, so 
that the conditions under which the individual trajectories were processed changed with the 
identity of the trajectory. We assumed that the different trajectories were emitted from the 
star successively, but ceased t o be emitted when th e collapse to the black hole occurred. 
This would be consistent with I Woosley et al.l (119941 ). who assumed successive emissions of 
the trajectories to generate what turned out to be a good representation of the Solar r-process 
abundances. 

For our r-process calculations, we used a network code ba sed on libnucnet, a libr ary of 
C codes for storing and managing nuclear reaction networks (IMeyer fc Adams! 120071 ). The 
nuclear and reaction data for the calculations were taken from the JINA reaclib database 



( Cvburt et al. 1120101) . We performed calculations for trajectories 24 - 40 in the lWoosley et al. 



( 1l994l ) hydrodynamics model. For each trajectory, reaction n etwork calculations w ere per- 
formed for T9 <2.5 using initial abundances derived from the IWoosley et al.l (119941 ) results. 
Our calculations were simplified by assuming an initial abundance of massive nuclei from 
a single nucleus heavy seed with a mass A equal to the average mass at Tg=2.5 and an 
atomic number derived fro m the average mass n umber, the Y e , and the neutron and alpha 
mass fractions at T 9 =2.5 in lWoosley et al.l ( 119941 ) . We justify this by noting that the heavy- 
nuclide distribution is typically sharply peaked near Tg=2.5. An adiabatic expansion was 
assumed for each trajectory in the nucleosynthesis, with entropy constant within a trajec 



tory b ut varying between trajectories, again consistent with the approach of IWoosley et al. 



( ]1994J ) for times at which the temperature T 9 <2.5. The materi al was taken to expa nd at 
constant velocity on a timescale consistent with that derived from lWoosley et al.l ( Il994l ). For 
each trajectory the calculation was continued until neutron capture reactions had ceased and 
the abundance distribution versus mass number had frozen out. Our representation of the 
full r-process abunda nces, shown as the dots in Figure [fl is not as good as that achieved by 
Woosley et al.l ( 119941 ). but our calculations do produce the basic features of the full r-process, 
namely the mass 130 and 195u peaks. 



The simulation of IWoosley et al.l ( 11994J ) that produced a good r-process representation 
began with trajectory number 24 and summed the nucleosynthesis yields from the remain- 
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ing 16 elements since all of these trajectories would have been ejected into the interstellar 
medium. We also began with trajectory 24, and performed a mass weighted sum of the 
nucleosynthesis from the trajectories up to some higher number trajectory to observe the 
resulting nucleosynthesis when the trajectories beyond our maximum trajectory were as- 
sumed to be consumed by t he collapse to a black hole. This produced a result similar to 



that of IWoosley et al.l ( 11994J ) when all trajectories from 24 to 40 were included. The results 
are shown in Figure [H There it can be seen that truncating the r-process at increasing tra- 
jectory number does terminate the r-process nucleosynthesis at increasingly higher nuclear 
mass. Note that although the curve representing trajectories 24 through 31 does reach the 
mass 195u peak, the abundance produced in that calculation is nearly two orders of magni- 
tude below that of the full r-process, which would render the higher mass nuclides difficult 
to observe. The abundances for that calculation would therefore appear observationally to 
terminate at a mass of about 140 u. 

In Figure |2] we compare several tr-process calculations with the derived elemental abun- 
dance pattern in the metal poor halo star HD 122563. This low-metallicity star ([Fe/H] = 
-2.7) is deficient in the heavy neutron-capture elements (Ba and heavier) relative to the light 
neutron-capture elements (Sr through Cd) when compared with the scaled Solar r-process 
pattern. The HD 122563 abundances are a very poor match to the scaled Solar s-process pat- 
tern. Its abundance pattern matches the scaled Solar r-process patt ern better up to an atomic 



mass of about 70 (mass o f ~174 u), but even this fit is unsatisfying (jSneden &: Parthasarathy 



19831 ; iHonda et al.ll2006l ). Stars like HD 122563 may be candidates for enrichment by the 
tr-process. Figure |5] demonstrates that the tr-process predictions, while far from a perfect 
match to the individual abundances, can reproduce the relative downward trend in abun- 
dance with increasing atomic number seen in some metal poor stars. More exploration of 
tr-process calculations is obviously needed, but the general trend is encouraging. 



4. Probability of Occurence of tr-Process Stars 

In principle, a test of the tr-process model would be provided by a large set of data 
for metal poor stars that spans the masses of the nuclides produced in the r-process from 
mass ~80 to the heaviest observable nuclides that the r-process synthesizes, usually thorium. 
Unfortunately this is challenging due to the difficulty of observing the higher mass rare-earth 
nuclides, starting at mass ~160. Thus all r-process events that terminate at masses between 
160 and lead would a ppear to terminate at mass 160, producing the effect observed by 



Roederer et all fl2010al ). 



Any standard initial mass function (IMF) reflects the fact that less-massive stars are 
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more numerous than more-massive stars. Assuming that stars from 8 to 40 solar masses 
may produce comparable amounts of r-process material, we can use the IMF to estimate the 
fra ction of s t ars w hose r-process may be truncated by collapse to a black hole. We adopt 
the lSalpeterl (Il954i ) IMF, for which dN/dm oc m -2,35 for massive stars, where m is the stellar 
mass in units of the solar mass, and N is the number of stars of a given mass per unit volume. 
The ratio of the number of stars that would be expected to collapse to black holes (25-40 
solar masses) to those expected to collapse to neutron stars (8-25 solar masses) is 0.13 in a 
well-sampled IMF. 

While a number of r-process rich stars (here taken to mean stars with [Eu/Fe] > +1.0) 
have been discovered and studied in detail in the last two decades, these stars constitute a 



relatively minor fraction o f all metal-poor stars. More unbiased samples (IMc William et al. 



19951 ; iBarklem et al.l 120051 ) find th at stars with [Eu/Fel > +1.0 comprise <10% of all stars 
with [Fe/H] < -2.0. Figure 11 of iRoederer et al.l fcOlOal ) suggests that stars with [Eu/Fe] 
< are candidates for enrichment by a truncated r-process. Difficulties in detecting Eu in 
metal-poor stars with [Eu/Fe] < make it even more challenging to estimate how numerous 
these stars are. Ba is more easily detected and may be used to represent Eu and other 
heavy ele ments in stars l a cking s-process enrichment. From the Ba abundances in the large 
survey of IBarklem et al.l (120051 ) we estimate a lower limit of ~55% of metal-poor stars as 
candidates for enrichment by a truncated r-process. This is significantly higher than the 13% 
derived above assuming a Salpeter IMF, yet Figure [2] demonstrates that a tr-process may be 
one way to produce the heavy nuclides observed in some metal-poor stars. The discrepancy 
in these fractions could indicate that additi o nal nucleosynthesis c h annels, such a s thos e 
pr oposed by, e.g., iQian fc Wasserburgl (120081 ) . iFarouqi et al.l (120091 ). IWanajo et al I ( 1201 ll ). 



or 



Nakamura et al.l (120111 ). together with the tr-process may contribute to the chemical 



enrichment of the early Galaxy. Given that massive stars with short lifetimes will dominate 
the chemical enrichment at early times, though, the tr-process may very well have been a 
major source of heavy nuclei at these epochs. 



5. Other r-Process Issues 

As a possible additional benefit of the tr-process, we note that the stars that were 
truncated at mass points 28 or less produced no nuclides in the mass 130 peak or beyond. 
This would have the effect of boosting the yields of the lighter r-process nuclides relative 
to the heavier ones, filling in the mass 110 to 120 u region. It has been known for some 
time that r-process calculations that produce a mass 195 u peak underproduce the nuclides 
in the 110-120 mass region. The r-process nuclides that were even lighter than those were 
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underproduced by an even larger factor, but this would also be consistent with distributed 
cutoffs of the r-process nucleosynthesis. Thus, again guided by the IMF, but along with the 
mass dependent cutoff times, t he tr-process might potentially solve th is troublesome aspect 
of current r-process analyses fcratz et al. 2007 ; Farouqi et al. 2009 ). Conclusions in this 
regard, however, must await more detailed analyses. 

An obvious test of the tr-process model could occur from renewed effort to observe the 
higher mass rare earth elements in very metal poor stars with low [Eu/Fe] ratios. If elements 
in that mass region can be observed, and the tr-process is correct, additional observations 
would be expected to map out the black hole collapse times over that region. While no 
doubt challenging, a relationship between observations and black hole collapse times would 
be of great interest. 
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Fig. 1. — r-Process calculations using the trajectories from IWoosley et al.l (1l994h . but 
summing the r esults from trajectory number 24 to some later trajectory, as was done in 
Woosley et al.l (J1994J). The successive curves include trajectory 24, 24 through 26, 24 through 
28, 24 through 30, 24 through 31, 24 through 32, and 24 through 40, where the emission times 
of the latest trajectory in each sum are as indicated. The sums to higher number trajectories 
make both the mass 130 and 195 r-process peaks, and t he sum through trajectory 40 come s 
closest to representing the Solar r-process abundances (IKappeler. Beer, fc Wisshak 1119891 ). 
shown as the dots. 
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Atomic Number 



Fig. 2. — C o mpar ison of tr-process predictions, using the formalism adapted from 
Woosley et al.l ( 11994 ). for four mass cuts. These are cuts below shell 30 (green), 31 (light 



blue), 32 (dark blue), and 40 (purple) with abun dances in the metal-poor star HD 122563. 
Observational data are from iHonda et all J2006h and Roederer et al. (2009, 2010b). The 
predictions are scaled to the La abundance (Z=57). 



